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Abstract 

The study of one or more baryons in lattice QCD is severely hindered by the exponential decay 
in time of the signal-to-noise ratio. The rate at which the signal-to-noise decreases is a function 
of the the pion mass. More precisely, it depends on the minimum allowed pion energy in the 
box, which, for periodic boundary conditions, is equal to its mass. We propose a set of boundary 
conditions, given by a "parity orbifold" construction, which eliminates the zero momentum pion 
modes, raising the minimum pion energy without altering the QCD ground state, and thereby 
improving the signal-to- noise ratio of (multi)-baryon correlation functions at long Euclidean times. 
We discuss variations of these "restless pions" boundary conditions and focus on their impact on 
the study of nuclear forces. 
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I. INTRODUCTION 



Lattice QCD studies of heavy systems are plagued by large statistical noise. The signal- 
to-noise ratio of a correlation function created by an operator with n quark-anti-quark pairs 
decreases with (Euclidean) time as g-C^-f where E is the mass of the state under 
consideration. This is a particularly nasty problem for the recent studies of nucleon-nucleon 

nn n 

[2j and hyperon-nucleon [3j forces with lattice QCD. The large statistical error renders 
the numerical information at large times useless. Compounded with this problem, at early 
times the correlators are contaminated by excited states and so there is only a very narrow 
range of time slices left containing useful information. In unquenched calculations 
the statistical noise allows for semi-quantitative results only, even after the computation of 
thousands of fermion propagators. These errors are also much larger than finite-volume 
and finite lattice spacing 5| effects in these observables. 

We propose here a scheme to alleviate this signal-to-noise problem. We start with the 
simple observation that the statistical noise is dominated by the energy of the lightest pion 
states. Periodic boundary conditions allow a pion zero mode and thus the lowest pion 
energy is equal to its mass. If one were to impose anti-periodic boundary conditions for all 
three pions, then the pion zero-modes are forbidden and the minimum energy is given by 
(assuming anti-periodic boundary conditions in all three spatial directions) 



7T J 



with L the size of each spatial direction. Thus it is clear that the signal can be improved 
by using these "restless pions" boundary conditions. There are other applications for which 
these restless pions are useful. In addition to the obvious benefit to spectroscopy studies 0] , 
the advantage of anti-periodic boundary conditions for the pion in the extraction of the 
K — > ixtx amplitude using the Lellouch-Luscher method ?| were pointed out in references 

0,3- 

In lattice calculations, one does not have direct control over the hadronic boundary con- 
ditions. What can be controlled at will are the boundary conditions of the quark and 
gluon fields. However, it is not obvious which modifications of the quarks and gluons at 
the boundary implies an anti-periodic boundary condition for all three of the pions; anti- 
periodic boundary conditions for the neutral pion have remained elusive. Obvious choices 
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get tantalizingly close to the desired restless pions but upon close inspection, have unde- 
sired consequences. For instance, twisted boundary conditions [k]] allow for continuous 
momentum transfer by providing hadrons a momentum kick at the boundary. But as we 
will explain in sec. [TTJ these twisted boundary conditions do not effect the signal-to-noise 



issue we are interested in. The G-parity boundary condition suggested in references 



breaks both the spatial subset of the hypercubic rotation invariance as well as chiral sym- 



metry ll|. In Ref. 9j, an isospin boundary condition was used, q(L) = r 3 g(0), but this 
leaves the neutral pion unaffected. This allows for an extraction of the A / = 3/2 K — > nit 
amplitude (in which the pions are in an I = 2 final state) but it does not serve our purpose 
of reducing the statistical noise for baryon calculations. Various "hybrid boundary condi- 
tions" have been employed, first in the numerical study of penta-qaurk states }12| |. In this 
first implementation, the u and d quarks were given anti-periodic boundary conditions while 
the s quark was given periodic boundary conditions, allowing for a definitive identification 
of bound vs. scattering states by forcing the ni^-system to have non-zero relative momen- 
tum (the scattering state) while leaving a possible ® + (uudds) resonant state unaffected (the 
bound state). In a second variation of the hybrid boundary condition, used in the numerical 

n n 

study of charmonium [13J and a possible tetra-quark state [14], an anti-periodic boundary 
condition is imposed upon the quarks while the anti-quarks are given periodic boundary 
conditions. This allows for the same identification of bound and scattering states for these 
systems as the first variant, however this second hybrid boundary condition violates charge 
conjugation invariance. Furthermore, neither variant of these hybrid boundary conditions 
will help with the signal to noise issue we want to address. 1 An axial twisted boundary con- 
dition, q(L) = 75^(0) (and similar choices) provide for anti-periodic pions but additionally 
make a ~ (qq) anti-periodic and alter the QCD pattern of symmetry breaking. 

We propose a novel approach to this problem making use of an orbifold boundary con- 



dition. Similar constructions have been employed in the context of the "chirality prob 



in extra-dimensional extensions of the Standard Model 1151, domain- wall fermions 
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em" 
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and the Schrodinger functional formalism 18j. Instead of relating the field values at the two 



ends of the box (z = and z = L), we impose periodic boundary conditions on an extended 
box — L < z < L. However, the fields at negative values of z are not independent, but are 



1 See section |n] for details. 
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FIG. 1: Log of signal-to-noise ratio of the two-nucleon correlator in the spin singlet channel as a 
function of Euclidean time (from the calculation described in [20j]). The pion mass is about 350 



MeV. The signal-to- noise estimate of eq. ([S]) was normalized to the lattice calculation at t = 11. 

determined from those with positive z. By appropriately choosing a relation between the 
quark and gluon fields in these two halves of the lattice (the orbifold condition), we can 
enforce a 7r(— z) = —tt(z) condition, eliminating the zero momentum mode for all pions, 
making them restless. 



II. SIGNAL-TO-NOISE RATIO ESTIMATES 



Here we review the argument estimating the statistical noise for lattice QCD correlation 
functions 19[. Consider first a nucleon correlator 

C(t) = (q(t)q(t)q(t) g(0)g(0)g(0)>, 

where, for clarity, we have suppressed the Dirac, flavor and color indices. At large times, 
C(t) is dominated by the intermediate state of lowest energy with the quantum numbers of 
the nucleon: 

C (t) ^ Ae~ M \ (1) 
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where M is the nucleon mass. In a Monte Carlo calculations, C(t) is estimated by an average 
over N gauge configurations 



C{t) = C(t) = ±Y, s A(t)s A (t)s A (t) 



N 

(s 3 A (t 



A 

^3 , 



where S A (t) is the quark propagator in each one of the gauge configurations, A. The variance 
in this estimate is given by 



^ = ^^2\S A (t)S A {t)S A (t) - C{t)\ 2 
A 

= (S\(t)S{ 3 (t))-\C(t)\ 2 . (3) 



For large times, C 2 (t) ~ e 2Mt , while the large time behavior of (S A (t)S A 3 (t)) can be found 
by noticing that 

(S A (t)S A 3 (t)) = (q 3 (t)Q 3 (t) g 3 (0)Q 3 (0)> 

^ Be~ M , (4) 

where Q is a fictitious quark with identical quantum numbers and properties of the q quarks. 2 
The long time behavior of the correlator in eq. (j3J) is then dominated by the intermediate 
state with the lowest energy with the quantum numbers of three qQ mesons. Since they 
have the same mass as the qq mesons, this lowest energy state is given by three times the 
pion mass. Thus, for sufficiently light pions, (S A (t) S A 3 (t)) decays at a rate smaller than 
C 2 (t). The signal-to-noise ratio of the nucleon correlator is then given by 

C ® AVN^- ~ y/Ne-W-t™** . (5) 



We show in fig. (FTT) the siernal-to-noise ratio in an actual lattice QCD calculation (details of 

nn 

the simulation can be found in Refs. [l|, l20|) as well as the estimate in eq. p]). 

The estimate in eq. (jSJ) is easily generalized for correlation functions of multi-baryons 
and baryons with strange quarks. In the case of two-nucleon correlators, for example, the 



2 This explains why the twisted and hybrid boundary conditions do not help the signal-to-noisc problem. 
The fictitious Q quarks have the same boundary conditions as the q quarks, and thus the qQ and Qq 
mesons have periodic boundary conditions and are allowed a zero-momentum mode. 
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signal-to-noise ratio is proportional to v / iVe~ < - 2M " 3m7r - ) *. Recent lattice studies of nuclear 
forces (and hyperon-nucleon interactions) were severely hindered by the fast decrease of 
the signal-to-noise ratio with time . The correlators at short times cannot be used for 
fitting purposes since it is contaminated by excited states 3 while at later times the statistical 
noise overwhelms the signal, leaving only a very narrow plateau from which the physics is 
extracted. This is in stark contrast to lattice calculations of tttt interactions 
two-meson systems j^j. 



21] and other 



III. PARITY ORBIFOLDS 



Let us now describe the basic idea of the orbifold construction in the case when only 
one dimension is orbifolded. Consider a lattice whose z coordinate belongs to the interval 
[0, L\. Extend it to [-L, L] and identify the points z = —L and z = L, effectively turning the 
interval [-L, L] into a circle. Let all fields, 4>(z), satisfy the periodic condition <f)(L) = <p(—L). 
Then identify the points z and —z by relating (p(z) to <p{— z), effectively transforming the 
circle into a line segment (including the boundary) as shown in fig. (T5]). In the simplest case, 
4>(z) = ±<p(—z). If the plus sign is chosen, <p(z) will be a linear combination of spatially 
symmetric wavefunctions, 

oo 

= £ 4° «>s(^). (6) 

n=0 

If, however, the minus sign is chosen then the 4>(z) will be a linear combination of anti- 
symmetric wavefunctions, 

^z) =f><."> sin (^), (7) 

n=l 

and consequently there is no zero mode for this field. The lowest momentum allowed is 
kmin = f with an energy of yj (n/L) 2 + m 2 . This upward shift in the minimum allowed 
energy value is the desired result. In order to eliminate the pions at rest we will require that 
7r(z) = — 7r(— z). The ways to achieve this by imposing orbifold conditions on the quark and 
gluon fields and the generalization to higher dimensions will be discussed next. 

3 In the single baryon sector, a significant improvement in the isolation of the ground state and excited 
states at early times has been achieved with the use of multiple operators combined with quark and gluon 
smearing 0] . The equivalent study for operators coupling to multi- nucleoli states has not been performed 
and is anticipated to be significantly more challenging and costly given the larger number of operators 
and quark contractions. 
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FIG. 2: Identification of z and — 2 points reduces the circle to a line segment. 
A. One-dimensional S 1 /^ parity orbifold 

In the simplest version of our proposal the orbifold trick is used in only one of the spatial 
directions. Consider QCD fields in the periodic box [0, L] x [0, L] x [— L, L] x [0, 0\ satisfying 
the "parity orbifolding" condition (the issues we discuss here belong to the infrared regime 
and we use a continuum notation) 

A»(t,x,y, z) = A^t, x,y,-z), for // ^ 3 

A 3 (t,x,y,z) = -A 3 (t,x,y,-z), 

q(t,x,y,z) = P g q(t,x,y,-z), 

q(t,x,y,z) = q(t,x,y,-z)V z , (8) 

where V z = 27573 is the ^-parity operator corresponding to a reversal of the z direction and 
we work in Euclidean space. 4 The ^-parity operator V z is obtained from the usual parity 
operator 70, corresponding to a simultaneous reversal of all three spatial axes, combined 
with a rotation by ir around the z-axis. The conditions in eq. (jHJ) relate the QCD fields in 
one side of the box to their parity conjugates in the opposite side. Notice that, since parity 
is a symmetry of the theory, the contribution to the action from the z < region is exactly 
the same as the z > region and the computational cost of using the extended box, [-L, L] 
is the same as that of the smaller box, [0,L]. The only effect of the orbifold condition is 
on the link connecting the z < and z > regions. In other words, it acts as a boundary 

4 We use the conventions 75 = 7^ = 1, 7^ = 7^. 
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condition at z = 0. In fact, consider the orbifolded action in the case of Wilson quarks 

S = k [£-1(73 - r)q\ - £1(73 + r)g_i] + a 4 (<Mi + g-ig-i) H 

= -2^(73 + r)P, gi + 2a 4 g igi + • • • , (9) 

where k is the hopping parameter, the index on the quark fields denotes the position in z 
(the remaining coordinates are implicit) and the dots denote the contributions from the two 
sides of the bulk, z > and z < (which are equal to each other). We see then that the 
orbifolded [-L, L] lattice is equivalent to a [0, L\ lattice with some extra terms residing at 
the boundary, as is the case with any boundary condition. 5 

Notice that we could have equally used the opposite z-parity operator, —V z , implementing 
a reversal of all three spatial axis followed by a rotation by — 7r about the z-axis. The 
difference between rotating by tt in the positive or in the negative direction amounts to a 2ir 
rotation which, for spin- 1/2 fermions, leads to a minus sign difference between V z and —V z . 
Physical observables, being quark bilinears, generally do not depend on this sign. As can 
be seen in eq. fl2]), however, the boundary terms are linear in V z and are able to distinguish 
between the choice in sign of V z . This shows that the orbifold condition breaks the z — > — z 
symmetry. 

The parity orbifold condition on the quark and gluon fields implies orbifold conditions 
for the hadronic fields. If we identify the pion field with the tt ~ 975T5 interpolating field, 
we see that it satisfies the desired 

7r(t, x, y, z) = -vr(t, x, y, -z) (10) 

orbifold condition. In fact, the same condition will follow if any other pion interpolating 
field is used like, for instance, n ~ qrqF^F^, since it depends only on the fact that the 
pion has negative intrinsic parity. In fact, all parity odd operators will satisfy a condition 
similar to eq. ( jTOl while the parity even operators will satisfy the analogue equation without 
the minus sign. In particular, the a field a ~ (qq) has a zero mode and the QCD pattern of 

5 Notice that, contrary to the continuum case, the boundary conditions in lattice field theory are already 
contained in the action. Different lattice action terms localized at the boundary imply different boundary 
conditions in the continuum and the relation between them is, in general, a complicated dynamical 
question. 
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symmetry breaking is not affected by the orbifolding procedure. The nucleon fields satisfy 



N(t,x,y,z) = -V z N(t,x,y,-z), 

N(t,x,y,z) = -N(t,x,y,-z)V z , (11) 

as can be seen using the interpolating field N ~ qQ T T2C^^q. In the non-relativistic domain, 
V z = 27573 reduces to cr 3 and the allowed modes for the nucleon are 

(A 

cos(^p) , n x ,n y ,n z = 0, 1, • • • 

W 

(12) 

(A 

sin(^fi) , n XJ n y = 0,1,- ■■ ,n z = 1,2,- •• . 



N(x, y, z) = e ^ x + l 2 L v < 




Notice that only spin up nucleons can be at rest. Consequently we can construct a spin 
triplet two-nucleon state, like the deuteron, with zero momentum but a spin singlet two- 
nucleon state will necessarily have a minimum momentum equal to ir/L. This asymmetry 
between spin up and down is a consequence of the breaking of the z — > — z symmetry 
discussed above. 

Unfortunately, the boundary term shown in eq. (Q is not 75-Hermitian and the fermion 
determinant is not positive definite. This makes simulations with dynamical quarks sat- 
isfying the parity orbifold condition impractical. However, this method is perfectly suited 
to implementation in the valence sector only, i.e. only on the propagators generated in 



the background of dynamical configurations. In refs. |23l . |24J |. it was argued that up to ex- 
ponentially suppressed corrections, for many channels of interest including baryon-baryon 
channels, different boundary conditions can be used in the valence and sea sectors of the 
theory, known as "partially twisted boundary conditions". Therefore, gauge configurations 
generated with sea quarks satisfying periodic boundary conditions can be used with valence 
quarks satisfying "parity orbifold" boundary conditions. Intuitively, the possibility of using 
different boundary conditions for sea and valence quarks follows from the observation that 
sea quarks can "notice" their different boundary conditions only if they propagate around 
the lattice. But, for observables without annihilation diagrams, the propagation of sea 
quarks around the lattice is suppressed by e~ mL , where m is the mass of the lightest hadron 
made of sea quarks or a mixture of valence and sea quarks. In our case, this is the pion 
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mass. This argument is better appreciated by looking at the graphs in fig. (J3J), which display 
examples of processes contributing to baryon-baryon scattering. Only diagrams containing 
a baryon-baryon intermediate state give rise to power law volume dependence (below the in- 
elastic threshold). These two intermediate baryons are made of valence quarks and therefore 
satisfy the orbifold boundary condition. We stress that the rate at which the signal-to-noise 
decreases is set by the valence nucleon and pion masses. 

The increase on the pion minimum energy has an additional benefit. With the exception 
of the relation between two-particle energy levels and the S-matrix, described by the Liischer 
formula, finite volume effects are suppressed by factors of e~ EltL . An increase on the value 
of E n is then clearly beneficial. This is specially important for the exponentially supressed 
correction to the Liischer formula where the suppression factor, formally of order e~ E,rL , can 
be sizable for realistic lattices and periodic pions with = m w {4]. These finite volume 
corrections can be estimated using an extension of chiral perturbation theory adapted to 
the case where valence and sea quarks obey different boundary conditions in the molds 



of 



23 



24 



25 



26 



27]. 



B. Three-dimensional parity T 3 /Z2 orbifold 

The method of the previous section can be generalized in order to remove the zero- 
momentum modes of the pions in all three directions, further improving the signal-to-noise 
ratio. The simplest generalization of eq. (jHJ) is 

A (t,r) = A (t,-r), 

Ai(t,r) = -A(t,r), for z= 1,2,3 

q(t,) = Vq{t,-r), 

q{t,r) = q(t,-r)V, (13) 

where V = 70 is the usual parity operator corresponding to the reversal of all three space 
directions. While the boundary conditions in eq. (jHJ) can be seen as a mirror placed at z — 0, 
the conditions in eq. ( fTBT) can be visualized as a pin hole located x = y = z = with a lattice 
[-L/2,L/2] x [-L/2,L/2] x [-L, L] x [0,/?]. Again, all three pions obey the odd orbifold 
condition ir(t, r) = —ir(t, — r) but now their minimum energy is v/3(f) 2 + m 2 n . Nucleons 
obey the same conditions as the quarks, N(t,r) = 7oiV(i, — r). Since in the non-relativistic 
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FIG. 3: Examples of two-nucleon graphs containing sea quarks. The left column shows the graphs 
at QCD level (dotted lines represent sea quarks) and the right column represents the same graphs 
at the low energy effective theory level. The graphs on the first row are proportional to e~ A ® CDL , 
the second and third are proportional to e~ m7lL . The last row shows a graph with a power law 
dependence on the volume. 

limit 70 reduces to 1, non-relativistic nucleons satisfy periodic boundary conditions and 
contain zero modes. This property is very convenient when extracting low-energy phase 
shifts on the lattice, as with the 3-D parity-orbifolding, there is no restriction on the spin- 
isospin channels one can study in the ground state and the standard Luscher formula relating 
energy levels to phase shifts is unchanged. As it will be exemplified below, the increase in 
the signal-to-noise ratio is dramatic. 
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(a) (b) 

FIG. 4: Left: Estimate of the signal-to-noise ratio with the S" 3 /^ orbifold condition 
g -(2M-3- v /( r ) 2 +m 2 )t (g \[(^ \[ ne ^ anc l with periodic boundary conditions e -( 2M -3m)t (daghed ii ne ) 

as a function of t. Right: Log plot of the signal-to-noise ratio with the T 3 /Z>2 orbifold condition 
g -(2M-3 A /3(f ) 2 +m 2 )t ^ go j ic [ j ine ~) anc j periodic boundary conditions e -(2M-3m)t (d as hed ii ne ) 

as a function of t. In both figures the pion mass is 350 MeV and the box size is L = 2.5 fm. 
IV. IMPACT ON LATTICE CALCULATIONS 



A. Nuclear force studies 



In order to provide an explicit example, we use the values of the parameters used in [l| 
to estimate the impact of the method advocated here in the expected rate with which the 
signal-to-noise ratio decreases with increasing time. We disregard the interaction energy 
between the hadrons and approximate the energy of the two-nucleon state by ~ 2M. The 
energy of the three-pions is approximated by ~ 3m n when periodic boundary conditions are 
used, 3^(f) 2 + ml if the S 1 /Z 2 orbifold is used and 3y/3(z) 2 + ml if the T 3 /Z 2 orbifold 
is used. The result is plotted in fig. (jl]). The inclusion of the interaction energy between 
the two nucleons would change the figure by very little. In fact, for pion masses above 350 



MeV the energy shifts found 



]J are of order of 10 — 20 MeV. It is expected, however, that 



28|, 



in a narrow band close to the physical value of m n the energy shift should be larger 
corresponding to the diverging scattering lengths, but still much smaller than the rest mass 
of the nucleons. Even the modest increase in the pion minimum energy found in the one- 
dimensional orbifolding has a potential significant impact by noise limited measurements. 
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In the case of the three-dimensional orbifolding that potential improvement is enormous 
(notice the log scale in the corresponding graph). 

B. Impact on K — ► tttt 
As pointed out in 

aa 

, the extraction of the K — > tttt amplitude with the Lellouch- 
Liischer method [j| can benefit from eliminating pion zero modes. The method to eliminate 
pions at rest discussed here can only be applied to the 1 = 2 channel. In the 1 = channel, 
the use of different boundary conditions in the valence and sea sectors alters the amplitude 
by factors that are not exponentially suppressed. Of course, a modified chiral perturbation 
theory taking into account the differences of the valence and sea sectors can still be used to 
relate the results of such a lattice calculation with the real world QCD amplitude. 

V. DISCUSSION 

We have introduced "restless pions" boundary conditions designed to reduce the rapid 
degradation of the signal-to-noise ratio which plagues studies of heavy systems with lattice 
QCD. We have shown how these boundary conditions can be implemented with a parity- 
orbifold construction in either one or three spatial dimensions. Unfortunately, the action 
at the boundary is not 75-Hermitian and so this particular construction is not suitable for 
the sea sector. However, this method is perfectly suited for implementation of the valence 
fermions. For non-scalar channels, the difference in sea and valence boundary conditions is 
felt only by exponentially small terms. The numerical cost of implementing these parity- 
orbifolded valence propagators is the same for propagators with ( ant i) -periodic boundary 
conditions, as the fields in each half of the bulk are not independent, and therefore the 
implementation is achieved with a special boundary condition on the non-doubled lattice. 
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